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Abstract. - It is shown that one can obtain quantitatively accurate values for the supercon- 
ducting critical temperature within a Hamiltonian framework. This is possible if one uses a 
renormalized Ifamiltonian that contains an attractive electron-electron interaction and renor- 
malized single particle energies. It can be obtained by similarity renormalization or using flow 
equations for Hamiltonians. We calculate the critical temperature as a function of the coupling 
using the standard BCS-theory. For small coupling we rederive the McMillan formula for Tc- We 
compare our results with Eliashberg theory and with experimental data from various materials. 
The theoretical results agree with the experimental data within 10%. Renormalization theory of 
Hamiltonians provides a promising way to investigate electron-phonon interactions in strongly 
correlated systems. 



Introduction. - For more than three decades Eliashberg theory [Q has been the standard 
theory to compute quantitative properties of superconductors like the critical temperature. For 
a review we refer to |^] . The basic advantage of Ehashberg theory compared to BCS-theory 
[Q is that it provides a framework to treat real phonon spectra and strong electron-phonon 
coupling. The critical temperature of a superconductor can be calculated with great accuracy 
[1). Eliashberg theory is based on a description in terms of Green's functions. It becomes 
exact in the limit where the electron mass is small compared to ion masses. 

In the past ten years, several new developments like the discovery of high-Tc materials 
and heavy fermions (for a recent theoretical review see 0|) stimulated research activities in 
strongly correlated systems |8| . Theories of strongly correlated systems are usually based on a 
Hamiltonian description, which is the natural basis to investigate properties of bound states. 
If one wants to investigate the effect of the electron-phonon interaction in strongly correlated 
systems, one has to set up a model Hamiltonian that contains phononic degrees of freedom. The 
main problem is that the standard Eliashberg theory cannot be applied to such models, because 
it cannot be formulated in a Hamiltonian framework. There are also other difficulties, for 
instance in heavy fermion systems, where the mass of the quasi-particles becomes comparable 
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to the ion masses, so that Ehashberg theory is no longer appHcable. Despite years of effort, our 
understanding of electron-phonon interactions in strongly correlated systems is rather poor. 
The major problem is that one needs a theoretical method based on a Hamiltonian description 
that allows an accurate treatment of the electron-phonon interaction. 

BCS-theory is based on a Hamiltonian approach, but it does not contain the electron- 
phonon interaction explicitely. Most of the properties of superconductors can be understood 
using BCS-theory. But one of the problems of BCS-theory is that the characteristic energy 
scale that determines the critical temperature Tc of the superconductor or the gap cannot 
be calculated. The famous BCS-formula = 1.130 exp(—l/Af(e^)F), where 9 is the Debye 
temperature and V is the strength of the interaction, has to be used to determine the parameter 
N{eF)V from Tc- Usually it is argued that the problems of BCS-theory are due to the fact that 
retardation effects in the interaction are neglected This viewpoint is possible if one has 
a theory in terms of Green's functions in mind like the Eliashberg-theory. In a Hamiltonian 
framework one usually works with instantaneous interactions. From a viewpoint based on a 
Hamiltonian the problem of BCS-theory is that the effective Hamiltonian has to be obtained 
from an initial Hamiltonian containing electrons and phonons. The effective electron-electron 
interaction Bardeen et al. had in mind was the phonon-induced interaction of Frohlich 
or Bardeen and Pines | [lO| . These interactions are constructed perturbatively. But regarding 
the energy scales in the problem it is clear that a perturbative approach must fail. The initial 
Hamiltonian is dominated by the electronic energies, which are typically of the order of a few 
eV. Further it contains phonons with an energy scale of the order lOmeV. Superconductivity 
arises due to a marginal relevant operator, the phonon-induced electron-electron interaction. 
The relevant energy scale in this problem is set by the critical temperature, which is typically 
at least an order of magnitude smaller than the phonon energies. Thus, if one wants to 
obtain accurate values for the critical temperature, one has to resolve an energy scale which is 
about five orders of magnitude smaller than the typical energy scale in the initial Hamiltonian. 
This means that one needs a renormalization procedure which allows to calculate an effective 
renormalized Hamiltonian from the initial Hamiltonian to high accuracy. 

Construction of the renormalized Hamiltonian. - A general renormalization procedure for 
Hamiltonians has not been available for longtime, although in his first paper on renormaliza- 
tion Wilson treated a Hamiltonian problem ||ll| . The only accurate renormalization procedure 
for Hamiltonians was numerical renormalization, which could be applied to fermionic single 
impurity problems. Three years ago, a new method for calculating effective Hamiltonians has 
been proposed by Wegner [^ . It uses continuous unitary transformations to calculate an 
effective Hamiltonian from a given initial one. This method has been called flow equations 
for Hamiltonians. It has been applied successfully to single impurity problem s [|l3| and to 
dissipative quantum systems |lj]. Independently of Wegner, Glazek and Wilson |15(| proposed 
a similar method, similarity renormalization. They use continuous unitary transformations 
to renormalize a given Hamiltonian. The aim of Wilson et al. is to treat light-front QCD 
[|6|. Recently Brisudova et al. calculated quarkonium spectra using this method ||l^. Both 
methods have a large range of possible applications in many particle physics. The aim of 
this work is to show that reliable quantitative results can be obtained with the framework of 
renormalization of Hamiltonians. 

Both methods, flow equations and similarity renormalization, have been applied to the 
electron-phonon problem ||l^, 19 1. Lenz and Wegner used the flow equations to calculate the 



effective electron-electron interaction. Their result differs significantly from the one obtained 
by Frohlich ^ or Bardeen and Pines . The interaction within a Cooper pair they obtained 
has no singularity and is attractive in the whole parameter space. In p9[ | I applied the similarity 
renormalization scheme to the electron-phonon problem. The effective interaction one obtains 
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is similar to the one calculated by Lcnz and Wegner. Furthermore, the single particle energies 
in the original Hamiltonian are renormalized. I calculated the critical temperature for a simple 
Einstein model. The aim of the present work is to show that Tc can be calculated for realistic 
phonon spectra using this method. The results agree very well with results from Eliashberg 
theory and with experimental data. 

Our starting point is the standard Hamiltonian for the electron-phonon problem 

H = Ho+Hi (1) 

with 

Ho = ^efc : clck : + ^Wg : 6^6, :, (2) 

k q 



Hi = ^{gk,qclck+qh\ + g*k,qc{^qCkhq). (3) 
k,q 

cj. and Cfc are the creation and annihilation operators for electrons. I have not included 
spin and band indices, but this can be done without difficulty. &J and hq are the creation 
and annihilation operators for phonons, here as well different acoustical and optical branches 
can be introduced. The colons denote normal ordering. In [ p9[ I showed how this problem 
can be treated using similarity renormalization. The basis of this approach is a continuous 
unitary transformation applied to the Hamiltonian. The transformation can be written in the 
form dH\/d\ = [ri\,H\]. A is a ultra-violett cutoff, and the generator rjx is chosen so that 
off-diagonal matrix elements vanish in H\ if the corresponding energy difference is larger than 
A. Details of this method in the present context are explained in a general description can 
be found in |15, During the transformation the electron-phonon coupling is eliminated 
successively and an effective electron-electron interaction is generated. The final result is an 
effective Hamiltonian that contains an electronic part of the form 

K\ ^^^k- clck ■■ ^'^'^'1 ■ A+qA'-^qCk'Ck :, (4) 

k k.k' ^q 

a part describing the phononic degrees of freedom and a weak electron-phonon coupling. 
The remaining electron-phonon coupling contains a small part of the initial electron-phonon 
coupling and other couplings involving two or more phonons. In jTot these couplings have 
been neglected since they are of higher order in the coupling constant gk,q- In principle it is 
possible to eliminate these couplings as well. This yields an additional contribution to the 
induced electron-electron interaction that is also of higher order. 



Gap and critical temperature. - In the following we will analyse the properties of the 
electronic subsystem described by the renormalized Hamiltonian (^) using BCS-theory. To do 
this we use the well known BCS gap equation 

^^=^^i^tanhf^ 

r2 , A2 V 2 



,r2 I A 2 
fct , ^ -I- 

It contains the interaction T4^_fe^g of two electrons forming a Cooper pair. In standard 
BCS-theory this interaction is often approximated by a constant in a small energy interval 
around the Fermi surface. The result of Frohlich for this interaction is attractive in 
a small region around the Fermi surface, it has a divergency due to a vanishing energy 
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denominator and becomes repulsive in the rest of the parameter space. In contrary, the 
result by Lenz and Wegner is attractive in the whole parameter space, it is of the form 
Vk,-k,q = l5fc,gP'^g/((efe — Cfc+ij)^ + ^g)- Similarity renormalization yields a similar result, 
Vk,~k,q — l5fc,gP/(|eI- ~ ffc+ql + Wq) • In the isotropic casc, and introducing the renormalized 
density of states N{(), the sum in (||) can be replaced by an integral. Furthermore, the 
electron-phonon coupling can be described by the standard Eliashberg function 2a'^F{uj) 
The gap equation can then be written in the form 



A(e) = / de' , ' ' tanh - sqrte'^ + A(e') 

The factor l/N^ep) is due to the definition of a'^F{uj), which contains a factor N{eF)- As a 
further approximation we replace the electronic density of states N{e) by a constant iVo. Since 
the electronic energies that enter in are renormalized, this has to be taken into account. 
As in Eliashberg theory PI, s imilarity renormalization yields a renormalization of the electron 
energies e^. — efe/(l + A) [|l9[ close to the Fermi surface. Therefore the renormalized density 
of states N{ep) at the Fermi surface can be replaced by N{eF) — No{l + A). A is the usual 
coupling strength in the theory of superconductivity, defined as A = (hj2a'^ F [uj) / to . In a 
first step we use (||) to determine Tc for weak coupling. If T is close to Tc, one can replace 
V'e'2 + A(e')2 by |e'|. We let e = and use that A(e') duj2a^F{uj)/ {e' + oj) is bounded, 
monotonic decreasing, and nonnegative. The second mean value theorem then yields 



A 



1 + A 



(ietanh(e/2Tc)/e = 1. (7) 



Keeping the leading logarithmic singularity this yields Tc oc cijexp(— (1 + A)/A). Cj can be 
obtained if one uses the second mean value theorem for tanh(e/2rc) and for A(e). Then one 
obtains Cj cx uJxog, LO\og = exp((lna;)), where the average (Inw) is taken with respect to the 
weight function a'^F{uj)/uj. Taking into account that the density of states N{e) is not constant, 
one can introduce an additional constant in the argument of the exponential function. In this 
way one obtains a Tc equation of the form 

Tc = c'cjiog exp(-ci^^). (8) 

It has the same form as the McMillan equation ||2^. For small coupling (A < 1) Tc/cjiog does 
not depend on the details of the phonon spectrum; this has also been shown in the framework 
of Eliashberg theory Q . 

In order to obtain results for Tc which can be compared with experimental data, one has 
to include the Coulomb repulsion. In principle this can be done form the very beginning. For 
our purpose it is sufficient to include the pseudo Coulomb potential |3[ || /x* by hand. With 
the assumption of a constant density of states the gap equation (^) can then be written as 

A(e) = I de' — ^ ^^"'^ tanh ( ^^e'^ + ^{e' f 



1 + A 



duj- (9) 



For small coupling one can again derive a McMillan-type formula. For small A the critical 
temperature does not depend on the details of a^F(a;), but only on wiog. To obtain accurate 
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Fig. 1. - Tc/tJiog plotted versus A for fi* = 0.1 and various types of spectra. The solid line is calculated 
for a lead-type spectrum, the long-dashed line for a mercury like spectrum, and the short dashed line 
for an Einstein spectrum. Also included is the McMillan curve (dot-dashed) and some experimental 
data. 



results, the gap equation has to be solved numerically for a given form of the phonon 
spectrum. In Fig. 1 we show numerical results for Tc in units of LU\og as a function of A for 
various types of spectra and ^* = 0.1. The solid curve is the results for a spectrum of lead. 
Here a'^F{uj) can be well approximated by a suitable sum of Lorentzians [Q. This has the 
advantage that the w-integral in (||) can be calculated analytically. The remaining integral can 
easily be calculated numerically to obtain Tc- The dotted curve is the result for an Einstein 
spectrum and has already been shown already in [l9|] . The dashed curve shows the result 
for a spectrum of mercury type. We have chosen these spectra since Allen and Dynes 
calculated Tc for the same spectra in the framework of Eliashberg theory. This allows a direct 
comparison with Eliashberg theory. The corresponding curves lie slightly above our results, 
the difference is about 5%. We have included in our plot some experimental data also shown 
in Q . A similar calculation can be done using the effective phonon-induced electron-electron 
interaction obtained by Lenz and Wegner |18 . The curves differ about 2% from the curves 



shown in Fig. 1. For the Einstein model, this has already been observed in ||19|. One notices 
that our results are systematically a few percent too small compared to the results from 
Eliashberg theory or to experimental data. The reason is probably that we neglected higher 
order terms in the effective Hamiltonian that describe an interaction between electrons and 
two phonons. Treating these terms in the same way as the electron-phonon interaction, one 
obtains an additional contribution to the effective electron-electron interaction, which is of 
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fourth order in the initial electron-phonon coupUng. This contribution is again attractive and 
will lead to a somewhat higher value of the critical temperature. 

Conclusions. - Our result shows that similarity renormalization or flow equations for 
Hamiltonians yield an effective Hamiltonian that contains the correct energy scale. The 
quantitative results are comparable to Eliashberg theory. But these methods work in a 
Hamiltonian framework, using continuous unitary transformations, and they do not rely on 
special properties like a small ratio of electron mass to ion masses. Therefore one may hope 
that the analytical treatment of electron-phonon interactions in strongly correlated systems 
is as well possible using these methods. 

From the viewpoint of Eliashberg theory or of a field theoretical approach based on a 
Lagrangian, where the phonons can be integrated out explicitely, a problem may still be that 
the effective interaction we calculated contains no retardation effects. But it should be pointed 
out that a quantity like or the energy gap is not a dynamical quantity, but a spectral property 
of the Hamiltonian. Dynamical correlation functions of some observable can be calculated as 
well using continuous unitary transformations. One has to take the transformation of the 
observable into account. This has been shown in a recent investigation of dissipative quantum 
systems in a Hamiltonian framework p^ . Using flow equations for Hamiltonians one can 
obtain accurate quantitative results for dynamical low temperature correlation functions as 
well. 

I gratefully acknowledge useful discussions with Franz Wegner, Stefan Kehrein, Stan Glazek, 
and Peter Lenz. 
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